Smearing of phase transition due to a surface effect or a bulk inhomogeneity in 

ferroelectric nanostructures 
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The boundary conditions, customarily used in the Landau-type approach to ferroelectric thin films 
and nanostructures, have to be modified to take into account that a surface of a ferroelectric (FE) 
is a defect of the "field" type. The surface (interface) field is coupled to a normal component 
of polarization and, as a result, the second order phase transitions are generally suppressed and 
anomalies in response are washed out. In FE films with a compositional (grading) or some other 
type of inhomogeneity, the transition into a monodomain state is suppressed, but a transition with 
formation of a domain structure may occur. 
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Theoretical studies of phase transitions in thin films 
and the corresponding size effects within the Landau the- 
ory [1,2] have been undertaken since 1950s. Recently the 
interest to these questions has risen dramatically in view 
of the applications of ferroelectric thin films [3] and a 
discovery of various ferroelectric nanostructures [4] . The 
boundary conditions for thin films were originally dis- 
cussed by Ginzburg and Landau (GL) in 1950 [5] and by 
Ginzburg and Pitaevskii in 1958 [6]. It was shown by 
GL that, if the properties of the boundary layer are the 
same as of the bulk, one arrives at the condition that 
the gradient of the order parameter vanishes at the sur- 
face, Vfj?7 = (in zero magnetic field, n is the normal 
to the surface). Starting from a microscopic theory de 
Gennes has shown that for a superconductor-metal inter- 
face with no current and magnetic field a more general 
boundary condition applies, Vfj?y + rj/6 = 0, where 8 is 
the characteristic length scale describing the proximity 
effect [7]. These conditions are very general and were 
obtained phcnomenologically by Kaganov and Omelyan- 
chouk for a surface of a ferromagnet [8] (cf. review in 
[9]). Kretschmer and Binder [10], using the same bound- 
ary conditions, have taken into account the depolarizing 
field, which is important when a ferroelectric polarization 
(or magnetization) is perpendicular to the surface. Later 
these boundary conditions have been used customarily in 
studies of phase transitions in ferroelectric films (see, e.g. 
[11])- 

It is obvious, however, that while the treatment [8,10] 
is appropriate for magnetics, it overlooks an important 
specific feature of ferroelectric phase transitions in thin 
films, wires, and other systems with boundaries. Indeed, 
there is an effective field at the surface (interface) of any 
material appearing because the surface breaks the sym- 
metry of the bulk. For instance, a part of this surface 
field might be due to a Coulomb dipole field (double 
layer), contributing to the work function [12,13]. This 
makes ferroelectric surfaces qualitatively different from 
the surfaces of magnetics. The effective field is coupled 



to the component of polarization perpendicular to the 
surface/interface and, as a result, the second order fer- 
roelectric phase transitions are typically smeared out, as 
we shall see below. 

We shall discuss, as an example, a paraelectric- 
ferroelectric phase transition in cubic perovskite thin 
films where a surface is perpendicular to one of the cu- 
bic axes. The surface or interface eliminates all the 
symmetry elements, which change a vector perpendic- 
ular to the surface and generates a local field conjugated 
to the polarization component perpendicular to the sur- 
face (P z in our case) [14]. To demonstrate the effect, 
we consider the state with P x , P y = 0, described for 
given potentials on electrodes [15,16] by the free energy 

F f = F LGD + fdV^-j: a e a Va, 
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where Vj_ = (d/dx,d/dy) is the gradient in a plane of 
the film, q a (ip a ) are the charges (electrostatic poten- 
tials) at the electrodes a = 1,2. Here Vj_ = for the 
monodomain state. We assume ideal electrodes with a 
vanishing Thomas-Fermi screening length. As discussed 
above, the surface produces an effective surface field w, 
and, generalizing Ref. [8], we have to add the surface 
energy to (1) to obtain the free energy of the film 



F = F f + 
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wP 



(2) 



where a corresponds to a "temperature" -like component 
of the surface energy. We obtain from Eqs. (1) and (2) 
after an integration by parts the correct boundary condi- 
tions for ferroelectrics 
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One can estimate that a ~ d a t, where d a t is the charac- 
teristic "atomic" length scale, on the order of the lattice 
constant. The electric field at the surface, first consid- 
ered many decades ago [12], is on the order of w/d a t ~ 
<b s /d at <~ IV/ A w 10 8 V/cm, where q$ s - 4cV is the 
typical workfunction for fcrroelectrics [3]. The surface 
bias field corresponds to a surface charge ~ 100/iC/cm 2 , 
which is on the order of an "atomic" polarization P at = 
q/d at ~ 200^C/cm 2 , so that w ~ P a td a t (we expect that 
the non-Coulomb contribution to w is of the same order 
of magnitude). The polarization P (z) is found from the 
equation of state for (1) and the Poisson equation, as- 
suming that there is no external charge, and neglecting 
for a moment the non-linear terms in polarization: 



d 2 P 



— {E + 4ttP) = 0, 

dz 
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j=E , 



(4) 
(5) 
(6) 



where E is the external electric field. We obtain from 
Eqs. (5) and (6) 



E = E - 4tt [P (z) - P] , 



(7) 



where the overbar means an average over the film, i.e. 
f = (1 /I) f_y 2 dzf(z). Substituting this into Eq. (4) and 
integrating over the film, we find 



Ap _g ( dPm _ dP(-l/2) \ = ^ 
I \ dz dz J ' 



(8) 



We write down the solution as a sum P = P + p( z ) of 
the homogeneous, P = (E Q + 4-kP) / (A + 4-rt) , and the 
inhomogeneous, p(z), term 

p = Ci exp [-\(z + 1/2)} + C 2 exp[-A(//2 - z)], (9) 
Cl(2) = (™1(2) - a 1(2 )P ) / ("1(2) + Aff) . (10) 



where A 
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g 



^ ~ d~ 



Since P = P + {d + 



C2)/\l, we obtain with the use of Eqs. (10), (8) 



A'P = E' + 



Wi 



w 2 



I \cti + \g «2 + A.g 



(11) 



where A' = A{1 - + \g£/l w A + 4tt£/AZ, E 1 = 
E Q (1 - S/M) , and f = + Xg) + a 2 /(a 2 + Xg) - 1. 

The phase transition in this case is smeared out, since 
generally the surface dipoles are asymmetric. In the sym- 
metric case, wi — —w 2 , a\ = a 2 , the phase transition 
persists, but the transition temperature of a transition 
into a monodomain state is shifted by the amount in- 
versely proportional to the thickness of the film, with the 
following estimate for displacive systems (cf. Ref. [10]): 



AT, = 



±<T at 1 
A V 



(12) 



The monodomain transition in the symmetric case occurs 
at A = —Xg^/l ~ —d at /l. This is close to a transition 
with the formation of domains [17]. Which transition 
actually occurs depends on materials parameters. 

The surface dipoles discussed above are a special case 
of polarization due to gradients of a scalar quantity (con- 
centration c, density, temperature, etc.) and they are 
accounted for by a term like 



f c = - 7 PVc, 



(13) 



in the free energy, where the coefficient 7 is estimated as 
7~ Patdat [18,19] (see also [20]). 

Consider now the case of a film with a compositional 
profile (grading) given by e.g. the concentration of one 
of the components of a ferroelectric alloy c = c(z). Such 
systems are currently a focus of research in ferroelectrics 
due to their unusual pyroelectric characteristics [21]. The 
equation of state of the graded ferroelectric film is 

d 2 P -de 
A{z)P + BP 3 - g^ - DV\P = E Q + 4tt(P - P) + 7 £. 

(14) 

Consider a special case of a step-wise concentration 
profile, i.e. c = C\ when < z < l\, and c — c 2 when 
— 1 2 < z < 0, and the boundary conditions are "neutral" 
(dP/dz = at z = l\, —l 2 ). The equation of state in this 
case is 

A r P + BP 3 -g^ = E Q + 4n(P-P), (15) 

for the both parts of the film r = 1,2. The boundary 
conditions at z = follow from the continuity of a dis- 
placement field E + 4irP and the equation of state (14). 
In displacive systems the electric field E ~ AP is al- 
ways much smaller than the polarization P, since \A\ -C 1 
[17]. Hence, with high accuracy oc (ci — c 2 )/Att <C 1, the 
boundary conditions are 



Pi =P2 

JdPi_dlV\ = _ h 
V dz dz J 



(16) 



at z = 0, with h = j(a - c 2 ) ~ P a td a t{ci ~ c 2 ). 

We have studied before [17] a similar situation but 
without the concentration gradient, i.e. for 7 = 0. Hav- 
ing assumed that the z— dependence of A is due to the 
concentration dependence of the Curie temperature with 
Td - T c2 - T at (ci - c 2 ) and dA/dT - T^ 1 (displacive 
systems) we have shown that for (ci — c 2 ) > d at /l the 
loss of stability of a paraphase occurs with a formation 
of a domain structure, and it takes place at A 2 > but 
A x < and \A X \ - (a - c 2 ) > d at /l [17]. For these 
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FIG. 1. Schematic of the polarization distribution in a fer- 
roelectric with a step-wise concentration profile. 

threshold values of concentration inhomogeneity the re- 
sult of Ref. [17] stands, whereas in the case of the bulk 
inhomogeneity and general nonsymmetric boundary con- 
ditions the results are different (see below). 

The main effect of the bias field is that now there is 
a polarization at all temperatures, and, therefore, the 
phase transition into a monodomain state is smeared out. 
However, a phase transition with formation of a domain 
structure is still possible. To see this, we need to in- 
vestigate a stability of a monodomain solution of (14). 
First, we need to find the average polarization P across 
the film. Integrating Eq. (15) over the film thickness, we 
obtain 

AP + BP 3 + SA r SP + B (sPW 2 + SP 3 ^ 

= E + w/l, (17) 
to = h + wx +w 2 - aiP(h) - a 2 P(-l 2 ), (18) 

where A = v\A\ + v 2 A 2 , 5A r = A r — A, vim = h(2)/l, 
I = h + h, Fig- 1- There are two possibilities: (i) near 
symmetric, \wi + w 2 \ 1^1 an d (h) asymmetric surfaces, 
l^i + w 2 \ 3> \h\. In the first case the monodomain transi- 
tion is smeared out by the presence of the gradient dipole 
field h, but the transition with domain formation is pos- 
sible. In the second, more general, case the monodomain 
transition is smeared out, and domains either form or 
they do not, depending on the concentration gradient 
and/or the thickness of the film. 

Importantly, the effective bias field w/l, conjugated to 
the average order parameter P, Eq. (17), is large. For a 
lOOOA thick film it would have the same effect as if there 
were an external electric field ~ 10 5 V/A, which, for com- 
parison, is only marginally smaller than the breakdown 
field in graded FE E b w 0.75 • 10 6 V/cm [21]. 

To investigate a stability loss of a paraphase, we as- 
sume that the linear approximation is valid and check 
later if the solution justifies the assumption. We es- 
timate P » w/ AI »j P at d at /(ci — c 2 )l, hence the first 



term in (17) is AP ~ w/l ~ P a td a t/l- The cubic term is 
BP 3 ~ P~ 2 (w/Alf ~ P at (d at /Al) 3 , since w ~ P at d at , 
and we can neglect it in comparison with the first linear 
term in (17) when A ~ (ci — c 2 ) > (d a t/l) 2 ^ 3 - The latter 
is the condition for the linearization of the equation of 
state (14), which takes the form 

(% + 4vr) SP - gd 2 SP/dz 2 = -w/l- SA r P, (19) 

with a solution 

r . , w/l + 5A r P , , . . 

SP r (z = L -^—+ Pr (z), 20 

A r + Air 

Pl = ae- XlZ + de- Xl(h - z \ < z < h, (21) 

P2 = be X2Z + C 2 e~ Xl(z+h] , -l 2 <z< 0. (22) 

We can replace Ai = A2 = A = y/4ir/g and obtain from 
the boundary conditions a — b = h/(2Xg) and 

XgCi + a±P(li) = w\, 
\gC 2 + a 2 P(-l 2 )=w 2 . (23) 

A reasonable approximation is P(h) k, C\ — 
wi/ (ai + Xg) , P(-h) ~ C 2 = w 2 / (a 2 + Xg) , Fig. 1. 
From first-principles calculations at the surfaces of 
BaTi0 3 and PbTi0 3 P - 10 _1 P ot (see, e.g. [22]). 

Consider the third and fourth terms on the left hand 
side of Eq. (1 7). Wit h the use of Eqs. (20)-_(22) we obtain 
the estimate 8A r 5P ~ (ci — c 2 )w/Anl AP, since there 
is an additional small factor (ci — C2)/47r. Both terms 
in SP r (z) give contributions to this estimate of the same 
order of magnitude. The term BPJW - P't 2 ^ (i^f ■ 
The condition that it is smaller than AP reads A ~ 
(ci — c 2 ) ^> (dat/Airl) 2 , and it is certainly obeyed when 
(ci — c 2 ) > (d a t/l) , which is the condition for the lin- 
earization, obtained above. The last term in (17) is very 
small if {dat/Aul) 2 <C 1, which is always the case. 

Now, we shall see if the domain formation is possible. 
Following the procedure of Ref. [17], we have to linearize 
(14) about the monodomain solution (inhomogeneous 
along z direction only) and look for its non-trivial solu- 
tions in the "soft" part of the film with A\ < in the form 
of the "polarization wave" , P(x,z) = P + 5P(z)+£(x, z), 
where £,(x, z) cx e lkx . We arrive at the same problem as in 
Ref. [17] but with a renormalized coefficient A\ — * Ai = 
A\ + 3BP 2 . The boundary conditions for £(x, z) are ex- 
actly the same as in Ref. [17], and A\ enters the condition 
for instability (A 2 > for the "hard" part does not). Do- 
main formation is possible when Ai ~ —d a t/l, in spite of 
a positive renormalization. The condition for this reads 
A {dat/1) 1 ^ 2 . This condition is stricter than the one 
for the linearity of the equation of state (17), meaning 
that our using of the linearized equation for P is justi- 
fied. There is also a range of concentration gradients, 
(d at /l) 2/3 < A ~ (a - c 2 ) < {d at /l) l/2 , when one can 
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linearize the equation of state for P, but domains do not 
form (A\ > 0). Finally, when the concentration gradient 
is even smaller, A < (d at /l) 2 ^ 3 , there is no domain for- 
mation and the equation for P is substantially nonlinear, 
AP < BP 3 . Therefore, the phase transition into mon- 
odomain state is smeared out, but a phase transition with 
the domain formation occurs when, in general, the con- 
centration gradient is large enough, c\ — c 2 > (dat/l) 1 ^ 2 , 
or if, for a given concentration gradient, the film exceeds 
some critical thickness, I > Id — d at j [c\ — c 2 ) 2 . If the 
system does split into domains in presence of the built-in 
surface bias field, the opposite domains will have different 
absolute values of polarization. 

In a special case of symmetric surfaces, when \w\ + 
w 2\ ^ the domains always form. Here we find that the 
net polarization is due mainly to concentration inhomo- 
geneity and is much smaller than in the general case con- 
sidered above, P « h/ Al ~ P a td a t/l- This means that the 
term BP 3 is on the order of P at (d at /l) 3 (since B ~ Pa~ 2 ), 
at the same time the term AP ~ (c\ — c<i)P a td a t/l, i- c - 
for A ~ (ci — c 2 ) > (dat/l) 2 the linear term indeed 
dominates in (17), AP > BP 3 . We obtain also that 
SP ~ h/Airl « AP/4tt ~ (ci - c 2 )P/4ir <C P. Because of 
this relation, all terms on the left hand side of Eq. (17) 
are indeed small in comparison with the first one, AP, 
and can be omitted. In the region of a stability loss 
with respect to domains A\ ~ —d at /l, therefore, 3BP 2 
~ (dat/l) 2 , i-e. the positive renormalization of A\ is 
very small, A\ = A\ + 3BP 2 rts A\ < and the system 
splits into domains. Therefore, for symmetric surfaces 
the presence or absence of the interfacial bias field at 
the boundary between two ferroelectric layers does not 
change our earlier prediction that practically any inho- 
mogeneity, however small, would lead to a domain for- 
mation [17]. 

One should make a reservation in case the bound- 
ary conditions correspond to a "surface ferroelectricity" 
(a < 0), then the monodomain transition can occur be- 
fore a domain structure forms. However, this is a some- 
what special case and, more importantly, effect of any 
real electrodes is rather opposite: it tends to suppress 
the ferroelectric transition into a monodomain state [23] . 
Therefore, it may be fairly difficult to observe the ef- 
fects of a "surface ferroelectricity" in the case of sponta- 
neous polarization normal to electrodes. Certainly, the 
real boundaries are never planar but rather rough. It 
seems likely that in real samples there are regions where 
the bias field is much smaller than would be in the case 
of planar boundaries. In these "weak" regions even small 
inhomogeneitics in the materials constants would lead to 
a formation of domains, just as in the case of a sample 
with the "neutral" boundary conditions. 

Note that even in a film with a step in a concentration 
profile the polarization is almost constant throughout the 
sample (with the exclusion of near-surface areas, Fig. 1). 



This is a result of a long-range depolarizing Coulomb 
field. It was neglected in a recent attempt to calculate 
the profile of polarization in the graded FE numerically, 
Ref. [24], and this led to erroneous conclusions. There- 
fore, those speculations do not apply to the observed be- 
havior, like a large apparent pyroelectric coefficient. In 
particular, the build-in bias voltage due to polarization 
inhomogeneity, that has been calculated in [24], should 
be exactly zero, 4ir J dz [P (z) - P] =0, Eq. (7). In fact, 
in graded samples there may be a bias voltage build-up 
due to a charge trapping, etc., and this could be related to 
the measured anomalous pyroelectric properties of these 
films, see e.g. a discussion in Ref. [25]. It is worth men- 
tioning that dipoles, introduced by interfaces, are likely 
to be important in ferroelectric superstructures, where 
they can affect an electric response of the structures [26] . 
We thank A. Kholkin for stimulating discussions. 
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